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Abstract 

We investigate the collider signatures of the multiple goldstini scenario in the framework 
of gauge mediation. This class of models is characterized by a visible sector (e.g. the 
MSSM or any extension) coupled by gauge interactions to more than one SUSY breaking 
sector. The spectrum consists of a light gravitino LSP, behaving as a goldstino, and a 
number of neutral fermions (the pseudo-goldstini) with a mass between that of the LSP 
and that of the lightest particle of the observable sector (LOSP). We consider the two 
situations where the LOSP is either a gaugino-like neutralino or a stau and we assume 
only one pseudo-goldstino of a mass of (9(100) GeV. The coupling of the LOSP to the 
pseudo-goldstino can be enhanced with respect to those of the gravitino giving rise to 
characteristic signatures. We show that the decay modes of the LOSP into a SM particle 
and a pseudo-goldstino can be significant. For both LOSP scenarios we analyze (pseudo)- 
goldstini production at colliders. Compared to standard gauge mediation the final state 
spectrum is softer and more structured. 
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1 Introduction 



In the coming years we will receive the final verdict from the LHC experiments as to whether 
low energy supersymmetry (SUSY) is present or not in nature. So far, both ATLAS and CMS 
experiments have not seen any signal and imposed strong bounds on the masses of colored 
superpartners (TJ |2j [3]. 

In this paper we will work in the general framework of gauge mediation, where the SUSY 
breaking effects are communicated to the observable sector by gauge interactions and gravity 
effects are subleading. A subclass of these models, namely those with a light neutralino and a 
low SUSY breaking scale, has already been constrained by the search for prompt photons plus 
missing energy [H El El E] ■ Given that the most straightforward searches have not yielded 
results, it is important at this stage to consider broadening the class of models to assure that 
we are not missing anything. 

In all generality, in a SUSY model based on gauge mediation the Lightest SUSY Particle 
(LSP) is the gravitino (we consider models preserving R-parity). The characteristic signatures 
of the model are then driven by the nature of the Lightest Observable-sector SUSY Parti- 
cle (LOSP) and by its decay properties. The two most common scenarios are those with a 
neutralino or a stau LOSP. 

Even within this general framework, it is important to keep exploring less conventional 
avenues to make sure that we get the most out of the current searches. One of the motivations 
for our work is to explore the theoretical aspects and experimental signatures of the goldstini 
scenario [91 [10] (see also [11] for previous work and [121 1131 EH EH ESI EH EE] for follow up 
work), in the context of gauge mediation as previously considered in [19]. 

This scenario arises when there are more than one SUSY breaking hidden sectors interacting 
with the observable sector via gauge interactions. Each sector gives rise to a goldstino, a linear 
combination of which is eaten by the gravitino while the remaining fields give rise to light, 
neutral, spin one-half particles, the pseudo-goldstini. For simplicity, in the following we are 
going to consider a set up where there are only two decoupled hidden sectors and hence only 
one pseudo-goldstino, but the analysis is trivially extended to the more general case and the 
qualitative features of the model are unchanged. 

It was shown in [19] that the pseudo-goldstino acquires a mass in the range 1-100 GeV from 
radiative corrections, while all supergravity effects can be neglected. If this mass is lower than 
that of the LOSP, the pseudo-goldstino becomes the Next-to-Lightest SUSY Particle (NLSP). 
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In this case one may face situations where the SUSY breaking scale is low, and thus the 
gravitino is an almost massless LSP, behaving as the massless true goldstino (denoted by G) up 
to a very good level of accuracy. Yet, the LOSP decays predominantly through other channels, 
namely those containing a massive pseudo-goldstino (denoted by G'). One of the consequences 
of this fact is that it softens the spectra of the final states, and we investigate those effects for 
the two most likely situations where the LOSP is either a gaugino-like neutralino or a stau. 

In the case of a neutralino LOSP (denoted by x) t ne most accessible signature arises when 
X is mostly a gaugino, promptly decaying into a photon and a (pseudo)-goldstino. The final 
state will contain photons and missing energy. A mostly gaugino x can al so decay into a Z 
boson and we will analyze these two situations in detail. It would also be interesting to consider 
the case where the x has a large higgsino component enhancing its decay into a higgsj^] For 
a single sector this was considered in [20], (see also the recent work [2T] where mono-photon 
signals arise directly for higgs decay). In the case of the stau LOSP (denoted by f), the lepton 
number conservation requires a stau to be created in pairs and the final state contains two taus 
and missing energy. 

Our analysis is independent from many detailed features of the spectrum of superpartners, 
in the spirit of simplified models. When forced to pick specific values for sparticle masses and 
mixing angles we will choose those used in common benchmark points [22] such as SPS8 (x 
LOSP) and SPS7 (f LOSP), in order to make a comparison with the standard gauge mediated 
models. Other LOSP scenarios, e.g. higgsino-like neutralino LOSP or colored LOSP suggested 
by general gauge mediation (23] , are also interesting to be applied within the goldstini framework 
and will be reported elsewhere. It should also be emphasized that although the points based on 
a minimal GMSB model have been constrained by the Tevatron [2U [25] and by the LHC [8], 
the multiple goldstini scenario eases the constraints due to the softer final spectrum. 

We will refrain from considering cosmological issues, just assuming a long enough lifetime for 
the pseudo-goldstino such that it safely escapes the detector (see e.g. [HJ [19] for considerations 
on the lifetime appropriate to our set up). 

The paper is organized as follows: Section 2 provides the theoretical background needed to 
describe the processes of interest. We discuss the spectrum and the couplings of the pseudo- 
goldstini to the Supersymmetric Standard Model (SSM). We mainly emphasize the differences 
from the well known couplings of the true goldstino with SSM fields. As already mentioned 
above, in the many-goldstini scenario, there are additional particles with masses in between 

1 Shorthand for Brout-Englert-Higgs scalar. 
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the LOSP and the LSP, the pseudo-goldstini. They come from other hidden sectors, but have 
couplings with SSM particles dictated by their nature of being goldstini in each hidden sector. 
The pseudo-goldstino has an interaction Lagrangian with the visible sector similar in structure 
to the one of the true goldstino, with however couplings that are no longer fixed in terms of 
the SSM parameters and the (overall) SUSY breaking scale F, but are essentially additional 
parameters. They can be computed in a specific set up, but can also be considered free as far 
as phenomenological models are concerned. They can easily be enhanced with respect to the 
goldstino couplings. Given that this is the most original feature of these models, we focus our 
interest on analyzing the dependence of the experimental signatures on these parameters and 
on the mass of the pseudo-goldstino for some appropriately chosen simplified spectrum for the 
observable sector. 

In section 3 we discuss the case of neutralino LOSP. We exemplify some of the features 
of these models by first considering processes that do not involve the intricacies of strong 
dynamics and can be treated analytically to a large extent. We first study the possible decays 
of the neutralino into (7G), ("fG'), and (ZG') pairs. (As mentioned before, in this paper 

we restrict our attention to a mostly gaugino x, and do not consider its decays into a higgs.) We 
then compute the production of neutralino and pseudo-goldstino in electron-positron collision 
(e + e~ — > xG'), which leads to the 7 +$ signature. The analytic result for the cross section has 
been used to test our FeynRules implementation but it is also of interest in future e + e~ linear 
colliders such as the ILC. We then simulate the neutralino-pair production, e + e~ — > xXj which 
gives the 77 +]fi signal, also of interest for the ILC. 

We compare how the total and differential cross sections change in magnitude and shape 
with respect to the standard results obtained in a single goldstino scenario. In particular, the 
total cross section of the single photon processes can be enhanced considerably for a pseudo- 
goldstino since the couplings of the latter are no longer strictly tied to the overall SUSY breaking 
scale. Furthermore, the mass and the coupling of the pseudo-goldstino alters the shape of the 
differential cross sections with respect to those of a nearly massless goldstino/gravitino, giving 
rise to a softer and more structured spectrum for the photons. 

We then proceed to the simulation for pp collisions of relevance to the LHC. We focus on the 
following exclusive processes without jets. First we discuss the photon(s) plus missing energy 
signals. As a more promising channel, the / + /~ + 77 signature is studied in detail, which 
is provided by pair production of sleptons, subsequently decaying to a lepton and a neutralino 
LOSP. While the total cross-sections are independent from the masses and couplings of the 
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pseudo-goldstino, the distributions of emitted photons and missing transverse energy do depend 
on these parameters and we illustrate this point by simulating various examples. 

In section 4 we discuss the case of stau LOSP. We first study the decay widths and branching 
ratios of the stau to (tG) and (tG 1 ). We then consider stau pair production in pp collisions, 
studying the r + r~ +$t signature. There is large SM background, but a satisfactory significance 
can be achieved by imposing some kinematical cuts. We find again that the pseudo-goldstino 
alters the shape of the differential cross sections with respect to the ordinary gauge mediated 
scenario. 

We briefly conclude in section 5. 

2 The two-sector model and goldstini couplings 

2.1 General formalism 

The purpose of this section is to derive the relevant couplings of the (pseudo)-goldstino. We 
consider a set up where there are two completely decoupled hidden sectors, each communicating 
to the SSM through gauge interactions. In other words, each sector is a model of gauge 
mediation on its own. (For a review of gauge mediation see Each sector will then have 

a goldstino, G\ and G2 respectively, coupling to the SSM particles in a way dictated by the 
universal low-energy effective action of the goldstino, applied to each sector. If F\ and F 2 are 
the respective SUSY breaking scales of each sector, and F = y Ff + F 2 , then the true and 
pseudo-goldstino are given by 

G=j(F 1 G 1 + F 3 G 2 ), (1) 

G'=j{-F 2 G 1 + F 1 G 2 ). (2) 

For definiteness, we will always assume F\ > F 2 . It has been shown in [19] that, while G can be 
considered massless to all effects, G' acquires radiatively a mass that is of the order of GeV if 
the scales in the two hidden sectors are of the same order. Moreover, when there is a hierarchy 
in the SUSY breaking scales F\ ^> F 2 , the mass is enhanced to 

Fi 

mc ~ — GeV. (3) 

F 2 

Since the exact expression for mg/ is very much model dependent, for phenomenological pur- 
poses it should be taken as a free parameter, recalling that it will reasonably fall in the 1-100 
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GeV range. We will always assume that G' is the NLSP, i.e. it is lighter than the LOSP 
otherwise its presence is virtually impossible to detect. 

The couplings of G and G' to the SSM particles can be derived by considering the couplings 
of G\ and G2, each as if it were the only source of SUSY breaking. 

The physical processes we are interested in involve only vertices with at most one (pseudo)- 
goldstino so we will not have to deal with the intricacies of the full (pseudo)-goldstino La- 
grangian. We start by considering the couplings to a generic observable sector whose chiral and 
vector multiplets we denote by $ 3 (0, ip a , F$) and W a 3 (X a ,A^,D) respectively (gauge and 
flavor indices suppressed). After discussing the general features in this language, we present 
the explicit couplings in the context of the SSM, which are slightly more involved due to the 
additional presence of Electro- Weak Symmetry Breaking (EWSB). 

Let us begin by considering the true goldstino G, whose coupling at the linear level is fully 
understood since the seminal work of the '70s |27J. It is well known that there are two equivalent 
ways of writing the linear coupling of G to the matter fields. One is the derivative coupling to 
the supercurrent 

1 

F' 

where, in the conventions of [28] 



C a = -(d lt G a J£ + h.c), (4) 



r = a u a^D v <f)* - ia^F^ + i— / =a v a p at i \F up + -^-=a^\D. (5) 

The other action is the non- derivative coupling obtained from the previous one by integrating 
by parts and using the equation of motion to obtain A a = J^| e o . m . 

C B = -hG a A a + h.c). (6) 
b 

Since the non-conservation of the supercurrent is entirely due to the presence of supersymmetry 
breaking soft terms, the expression for A a must be a function of the latter. 
Let us consider the most general soft SUSY breaking terms, namely: 

£ soft = --m A AA- -m* x XX- U{<f>,<f>*), (7) 

where U is a (at most cubic) gauge invariant function of the scalars, containing for instance 
the sfermion soft masses and the B[i term and m\ are the Majorana masses for the gauginos. 
The divergence of the supercurrent is now 

«V£ = A Q = ^=a^\ a F^ - i^KD - i> a ^ (8) 
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from which the coupling follows using ([6]). This expression is valid regardless of whether the 
gauge symmetry is spontaneously broken or not (see for instance [29] for a complete treatment 
in the MSSM). The form of ([6]) can also easily be derived from a superspace formulation of the 
broken SUSY theory, in terms of the goldstino superfield X. 

The derivative and non- derivative actions are completely equivalent, of course, but some 
issues are easier to investigate in one formalism than in the other. 

It turns out that it is the non-derivative action that is more directly generalized to the 
pseudo-goldstino case. This can be seen in various ways but perhaps the easiest argument 
comes by looking at the leading high energy behavior of the 2 — > 2 scattering amplitudes 
involving a goldstino. The action @ contains terms of dimension six and one would expect 
the tree-level unpolarized squared amplitudes to scale like s 2 where a/s is the center of mass 
energy. This is not what happens however since SUSY ensures that the leading order behavior 
cancels between the different contributions. This must be so since (jl]) is equivalent to ([6]) which 
contains terms of dimension at most five and yields a scaling of order s. 

The same scaling must occur for the pseudo-goldstino since, after all, at tree level it is 
a linear combination of two decoupled goldstini, but now the relative coefficients between the 
various terms in the action are no longer fixed by SUSY. Using the non- derivative action ensures 
that the high energy behavior is preserved. One can of course use the equations of motion 
"backwards" and rewrite the non-derivative action in terms of the same type of terms that 
appear in the derivative action but with different relative coefficients. In the process however 
one also picks up dimension six contact terms schematically like GXipip that once again cancel 
the s 2 behavior. Thus it is clearly more convenient to work with ([6]). 

In order to derive the couplings of the pseudo-goldstino, let us start at tree level with the 
two fields Gh {h = 1, 2) that are the goldstini of the respective hidden sectors contained in the 
superfields Xh- For each goldstino-gaugino-gauge boson vertex, we have 

1 J d ^X h W> 3 ^-WG^ (9) 

For the goldstino-fermion-sfermion vertex, we have 

9 9 

/TTv 777 
#g_m Xh Xl<M>1 D -^CW. (10) 

The terms ([9]) and (jTOl) correspond to the first and last terms in (JSJ) inserted into ([6]), the last 
one with U = m|0*0. We ignore further vertices that do not contribute to the processes of 
interest but a full treatment can be found in [30J. 



6 



Rotating to the G, G' basis, we obtain the couplings 



Xa^GF^ + Kx—^-Xa^G'F^ (11) 



and 

2 2 

mi mi , 

-j=rGip(j)* + K^GW, (12) 

where mA = m^i) + m A ( 2 ) and m| = m^^ + m^ 2 ^ , and the factors K\ and If^ are the ratios 
between the coupling of the pseudo-goldstino to the one of the true goldstino. Their expressions 
are 

(13) 
(14) 

where in general there will be different coefficients for each gauge group and matter multiplet. 
In some specific models, there can be relations or bounds between K\ and K^. However, in the 
phenomenological set up, one assumes that these two parameters are free. In any case, below 
we will use couplings of G' only to a very limited set of (s)particles. The interesting case will 
be when K\,K<$> 3> 1. Note however that both cannot be larger than F1/F2, hence they are 
limited to be somewhat smaller than 10 3 . (If F\jFi is too large, the above picture is no longer 
valid because the SUSY breaking vacuum of the second sector is destabilized [TP].) 



K x = - 


m A 


F 2 
Ft 


+ 


m A 


Ft 

iV 


K^ = - 


m 5(U 


F 2 


+ 


m m 


Ft 


ml 


Fx 


m 2 A 


F 2 ' 



2.2 Pseudo-goldstino couplings in the SSM 

We now specialize to the case where the observable sector is a SUSY extension of the SM. 
A couple of technical issues arise since one must also rotate to the physical bases of mass 
eigenstates of the observable sector after EWSB. More generally the (pseudo)-goldstini will 
also mix with the particles in the observable sector carrying the same quantum numbers but it 
is possible to treat this mixing to first order in the SUSY breaking parameter. 

Let us begin with the coupling of the neutralino to the G/G' and the photon. This involves 
an element of the neutralino mixing matrix, essentially proportional to how much the neutralino 
is the would-be photino. More precisely for the goldstino there is a factor (see e.g. [28J) 

a 7 = N*! cos 9 W + N* 2 sin 9 W , (15) 
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where N* x and N£ 2 are the mixing angles between the lightest neutralino and the Bino and 
Wino respectively. For the pseudo-goldstino, we can multiply the above by a factor K 1 that 
can be larger than one as discussed above. The relevant part of the Lagrangian is thus 

£ xG7 + £ xG , 7 = ^ X a^GF, v + ^^ X a^G%„. (16) 

The coupling of the goldstino and the pseudo-goldstino to the neutralino and the Z boson 
is slightly more subtle. Standard computations, such as in [31], use the goldstino couplings in 
derivative form. After EWSB, there are two such couplings: 

Ca = i^^x^a^GF^ - a ^^ xa ^d,GZ l/ + h.c, (17) 

where 

a Zr = —N* t sin 9 W + N* 2 cos 9w, (18) 

a ZL = N* 3 cos (3 - N* lA sin 0. (19) 

N± 3 and N^ 4 are the higgsino components of the lightest neutralino, and tan (3 is the ratio 
of the vacuum expectation values of the two higgs doublet^. The notation follows from the 
notable fact that the first and second term respectively couple only transverse and longitudinal 
components of the Z to a massless goldstino. 

Integrating by parts and using the equations of motion, one obtains the following terms in 
the non- derivative Lagrangian: 

_ az T m x + az L m z m z (a ZT m z + a Zh m x ) 
L xGZ 2y/2F TIf 

It might seem at first that the second term in the Lagrangian above cannot be reproduced 
using ([8]). This term appears because after EWSB there are off-diagonal mass terms involving 
the goldstino and the neutralinos. The mass eigenstates have then to be shifted by 0(1/F) 
terms mixing the goldstino with the neutralinos. Eventually the term above arises from the 
neutralino gauge couplings, and its precise form follows from using the equations leading to the 
EWSB vacuum (for a discussion, also including the pseudo-goldstino, see [T6]). 

The pseudo-goldstino couplings to the Z boson will be of the form above, with however 
some model dependent factors in front of each term. Since there are two independent terms in 



2 Although we do not consider the decay into higgs, we will nonetheless treat the decay into Z exactly by 
including the contribution of the longitudinal component. 
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( 120]) . there will be two independent coefficients and we choose to associate them directly to the 
rescaling of the az T and a Zjj coefficients^ 

r K ZT a ZT m x + K ZL a ZL m z 

1~> X G'Z = r- a r\u> 



2^2F 

1 „ _l_ K „ n „ m \ 

X^G'Z^ + h.c. (21) 



.m z (Kz T a ZT mz + K ZL a ZL m x 

1 72> A " ^" 



We now move on to the couplings of the stau. Even in this case the general formula 
( II 2 p requires some well known modifications after EWSB. Namely, the presence of a SUSY 
contribution to the tau mass makes the coefficient in the action depend on the difference of 
the masses in the multiplet. The large soft off-diagonal corrections to the stau mass matrix 
require rotating from the gauge eigenbasis f L , f R to the mass eigenbasis f , f'. For the couplings 
to the first two families these effects can be neglected. We focus only on the LOSP f = 
cos 9f tl + sin tr and write the Lagrangian as 



2 2 
mi — m 



Cf GT + Ct&t = — ^ — - (cos 9 f Gr L + sin 6 f G^t r ) t* 
2 2 

+ m -~ m - (R TL cos B f G'r L + K TR sin 6 f G'W R ) f* + h.c. , (22) 

where we have set Pl/rt = t L / R . Once again, we treat the coefficients K TL and K TR as free 
parameters. 



3 The case of the neutralino LOSP 

We begin our analysis by considering the case where the neutralino % is the LOSP. For this we 
use the Lagrangians (TTB|) . (f2"UI) and (|2"TI) . 

3.1 Decays of the neutralino 

The partial widths for neutralino decay into a photon and a (pseudo)-goldstino are: 

alml 

r( *^ G) = T^ (23 > 



3 One could have also chosen to rescale each Lorentz invariant term in (|20|) independently but this would 
obscure the comparison in the simulation. 
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We could also write F = y/Sm^^Mp in the denominators above, with the mass of the gravitino 
(i.e. the true goldstino) m 3 / 2 and the reduced Planck mass M p = 2.43 x 10 18 GeV. 

For a rather massive pseudo-goldstino, the factor between parenthesis can be significantly 
smaller than 1, though it is always of 0(1) barring any fine tuning of against m x . Hence the 
branching ratios can be of the same order if Kj = 0(1), or we can have a neutralino decaying 
almost exclusively to the pseudo-goldstino if K-y ^> 1. For instance we can have ~ 100 and 
then the decay rate of the neutralino is going to be 10 4 times larger with respect to a single 
sector scenario. 

In order for the neutralino to decay inside the detector, its total width cannot be too small. 
A rough order of magnitude of the bound is r to t > 10~ 16 GeV. For a single sector scenario 
and m x ~ 200 GeV, it would translate to \/F < 10 3 TeV (see e.g. [32]), but in our case the 
constraint is more flexible due to the presence of the factor. 

We now list the partial widths for the decay of the neutralino into a Z boson and a goldstino 
or a pseudo-goldstino. For the goldstino, using either f JT7|) or (120]) . one obtains the classic 
result [31] 

(2a| +4 )ml ( ml\ 4 
V{x^ZG)= K \ n ™ X U-^f • (25) 



32vrF 2 V m l 



For the pseudo-goldstino, the decay rate is given by 

r( * - ZG) = teL-J m ^ (26 > 

where = /3(^, -^-) with 

(3{a, b) = (1 + a 2 + b 2 - 2a - 2b - 2ab) l/2 (27) 



is the usual phase space factor and the spin summed and averaged amplitude squared is 
1 

2F2 



|-M| 2 = {{m 2 x - m 2 z )\2K 2 + K 2 ) + Qm G ,m z (m 2 x - m 2 z ) 2 K T K L 



+ m 2 G ,(m 2 x - m|) [(-4mJ - m%)K% + (-2mJ + m 2 z )K 2 L - 6m x m z K T K L ] 
- Qm G ,m z \m x m z {Kl + K\) + (m\ + m 2 z )K T K L ] 

+m G , [(2m 2 x + m 2 z )K 2 + (m 2 , + 2m 2 z )K 2 L + 6m x m z K T K L ] } , (28) 

where K T = K Zt cl Zt and K L = K Zh a ZL . The decay channel is open only when m G i < m x — m z . 

We note that, upon setting = K Zt = K Zl = 1, the decay rates to a pseudo-goldstino 
slightly differ from those to a massive spin-3/2 gravitino, as detailed for instance in 
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Figure 1: Branching ratios for the lightest neutralino decay into a (pseudo)-goldstino and a photon 
or Z at the SPS8 benchmark point. On the left the coupling of the pseudo-goldstino to the photon is 
gradually enhanced, while on the right the case of enhanced coupling of the pseudo-goldstino to the 
Z boson is shown. 



In Fig. [T] we plot the branching ratios for the lightest neutralino decay for varying pseudo- 
goldstino mass and for various values of K y , Kz T and Kz L - The remaining relevant SUSY 
parameters are taken at the SPS8 benchmark point as calculated by SOFTSUSY [31]. Namely, 
we have m G (= m 3/2 ) = 4.74 eV, m x = 139.2 GeV, N u = 0.99, N 12 = -0.031, N 13 = 0.124, 
iVi4 = —0.048 and tan/3 = 14.5. In general, the partial widths for decay into Z bosons are very 
much suppressed with respect to the ones for decay into photons, except when K Zt and/or K Zl 
are the only large factors. In the following we will not assume this. Also, the partial width for 
decay into higgses is negligible at the SPS8 point. See [16] for a different set up where this is 
not the case. 

It is obvious that the neutralino total width will always be exceedingly small, e.g. r tot ~ 
10~ 12 GeV for = Kz T = Kz L = 1, compared to m x , so that we can safely place ourselves 
in the narrow width approximation (NWA) in all processes of interest with an intermediate 
neutralino. In the next sub-sections we will discuss observable signatures involving the produc- 
tion and decay of neutralinos. Thus, the (differential) cross sections will be proportional to the 
square of the amplitudes for production of neutralinos, and to the branching ratios for their 
decay. 
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3.2 Goldstini and single photon production in e + e collisions 



We now perform an analytic computation with the purpose of highlighting the differences with 
respect to the single sector case and the role played by the extra parameters K\ and K^, 
characterizing the pseudo-goldstino couplings. For simplicity, in this subsection only, we stick 



to the case where the neutralino is a pure photino, i.e. a 7 = 1 and az T = az L 



0. 



There are three kinds of diagrams contributing to this process, as in Fig. [2j In the s- 
channel, the intermediate particle is a photon, since the neutralino is pure photino. Note also 
that there is another s-channel diagram with the outgoing arrows reversed. (We are using the 
two-component notation of [28].) In the t- and w-channels, the intermediate particle is either of 
the right- and left-handed selectrons. We assume here that the factor is the same for both 
selectrons. 

The couplings of the (pseudo)-goldstino have been reviewed in section [2] and for a pure 
photino neutralino they are (see ffTTj) and f|T2l) ) 

2 

Tfl Tf]~ 

C G , D K\- 7=7-7^ Xa^G'F^ + K*-^ G'W, (29) 



'2V6M p m 3 / 2 



V3M p m 3 / 2 



where 7713/2 = F/y/3M p is the gravitino mass, mg is the right and left selectron mass (we 
have neglected the mass of the electron with respect to the mass of the selectrons) and finally 



rri\. The other couplings needed for the computation can be gathered from e.g. [28] . 



Note that in the amplitude squared all diagrams actually interfere due to the mass of the 
pseudo-goldstino. The computation is thus more involved than the one for a massless goldstino 
J . The differential cross section for production of a photino and a pseudo-goldstino reads 



da P 



-\M\ 2 d$ 2 



(30) 




G/G' 




1>) 



G/G' 




Figure 2: All diagrams contributing to e + e — > xG/G': a) s-channel, b) i-channel and c) ti-channel. 
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with the spin summed and averaged amplitude squared 
3 2 r Klml 



e 2 ( Kim 2 

\ M \ 2 = qm-2 2 { i 2tu ~ + u ) + Zm G ,m x s + m 2 G ,(2m 2 x - t - u)) 

6M p m 3 / 2 I s 

j,2 4 f (* ~ ~ ™Z) , (u - m 2 G ,)(u - m 2 ) 2m G ,m x s 
+ K.m, oxo + ; 5^ — — + 



(t — ml) 2 {u — m 2 -) 2 (t — m 2 ~){u — m§) 

'm x {t-m 2 G ,)+m G ,{t-m 2 ) m x {u - m 2 G ,) + m G ,(u - m 2 ) 
- 2K x K^m x m~ — — - 2 ^ + — — 2 ±- 



m| u — m\ 



(31) 



where e is the electromagnetic coupling constant. We refer to e.g. [37] concerning kinematics. 
One can check that in the limit m G ' — > and K\ = = 1 this reproduces the result reported 
in [3Sj. Plugging \A4\ 2 in (13D1) and integrating over cos 9 we get the total cross section. 

It is interesting to note that the above amplitude is slightly different depending on the 
relative sign of the respective (real) Majorana masses of the neutralino LOSP and the pseudo- 
goldstino. This relative sign is model dependent. 

It should be stressed that the cross section scales with K\^jm\i 2 . In the spirit of [SU ESI 
[39j HOj ST] , LEP bounds on such cross sections can be translated into upper bounds on K/m 3 /2, 
or alternatively on if at a given value of m%/2- (Here K is for simplicity a common value for K\ 
and K^.) Roughly, we get K < 10 4_5 (m 3 / 2 /eV), which allows us some elbow room. Note that 
in the case of the single true goldstino production, i.e. with K\ = = 1, the cross section 
is very small, unless the gravitino mass is of order 10~ 5 — 10 -4 eV [36]. In the case of the 
pseudo-goldstino, instead, the cross section can be enhanced by the couplings K\ and while 
keeping the gravitino mass to standard values for gauge mediation scenarios, i.e. m 3 / 2 ~ eV. 

To exemplify the physics of this process, and in particular the dependence on the pseudo- 
goldstino parameters, in Fig. [3] we plot the total cross section as a function of the pseudo- 
goldstino mass for some values of the parameters K\ and K^. Here we consider a photino 
LOSP and take the masses as m 3 / 2 = 10 -9 GeV, m x = 140 GeV and m iR = m iL = 400 GeV. 

There is a destructive interference between the diagrams and thus the cross section for 
large turns out to be greater than the cross-section when both K\ and are large. We 
notice that rather large values of K\ and are required to obtain the cross section around 
O(10 2 ~ 3 ) fb with the eV order gravitino mass, while such large values are not favored by the 
stability of the SUSY breaking vacuum as mentioned before. 

One can also easily see that in these cases the emitted photons can be significantly softer 
than in usual gauge mediation scenario, since the pseudo-goldstino has a non negligible mass. 
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Figure 3: Total cross section of e + e~ — > \G' at y/s = 500 GeV (black lines) and 1 TeV (red lines) as 
a function of the pseudo-goldstino mass, for various values of K\ and K^. 

Moreover, similar to the discussions in jl2], the photon energy distribution can tell us about 
the masses of the neutralino and pseudo-goldstino as we will explicitly see in the next section. 
On the other hand, different K factors would give different photon angular distributions. 

All the results presented here can be obtained numerically running MadGraph 5 |43j simu- 
lations adapted to the (pseudo)-goldstino scenario (building on |4"4"]). having implemented the 
model using FeynRules [45j HQ |47] . This provides also a non trivial test of our software package, 
which we will use afterwards to simulate pp collisions. 

3.3 Goldstini and di-photon production in e + e collisions 

The production of two neutralinos will lead to a di-photon plus missing energy signature, which 
is evidence for the processes 

e+e" -> 77GG, e + e~ ->• 77GG", e+e" -> 77CG'. (32) 

The total cross section (<t^__j. ~ 177 fb at y/s = 500 GeV at SPS8) is similar to the single 
sector case, since the couplings that can be enhanced in the pseudo-goldstino scenario only 
appear in the decay of the neutralinos. However the photon spectrum, and in particular the 
edges of the energy distribution, is sensitive to the mass of G', both if the branching ratios are 
comparable or if the decay to G' is favoured. 
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In Fig. HJ the distributions of the leading photon energy (left) and of the missing invariant 
mass (right) for e + e~ — > XX — > 77 at = 500 GeV are shown. To obtain both plots, we 
applied a cut on the energy and the rapidity of the photons, > 15 GeV and \tj^\ < 2, as the 
minimal cuts for the detection of photons. In addition, we imposed the invisible invariant mass 
cut M; nv > 100 GeV to remove the SM (Z — > uu)jj background. The remaining background 
comes from the t-channel P^-exchange process, and this can be reduced by using the polarized 
e ± beams. 

Besides the reference point m G i = with K 1 = 1 (for which we have essentially two 
indistinguishable copies of a light goldstino), we take two different pseudo-goldstino masses, 
85 GeV and 125 GeV, with different couplings as in the following table: 



la. 


m G , 


= 85 GeV 


with K 1 


= 1 


[B( X lG') 


~ 0.2] 


lb. 


m G > 


= 85 GeV 


With Kry 


= 2 


[B( X lG') 


~ 0.5] 


lc. 


m G , 


= 85 GeV 


With Kry 


= 10 


[B(x lG') 


~ 1] 


2a. 


m G i 


= 125 GeV 


with K 1 


= 10 


[B( X lG>) 


~ 0.5] 


2b. 


m G > 


= 125 GeV 


with K 1 


= 100 


[B( X lG>) 


~ 1] 



Here, we keep Kz T = Kz L = 1, and hence the decay modes into Z are negligible; see also Fig. [T] 
for the branching ratios. We have also run the simulations for m G i = 10 GeV with = 1,2 
and 10, for which the branching ratio is respectively 0.5, 0.8 and 1, but we found distributions 
essentially overlapping with the massless one. 

The edges of the energy distributions allow to determine both the mass of the neutralino 
LOSP and of the pseudo-goldstino. A simple generalization of the massless goldstino case 
discussed in [3T] gives the following expression for the minimal and maximal energy of each 
emitted photon: 



In this case, it is difficult to determine the minimal edges due to the detector cut. On the other 
hand, the E™ x is rather clear although the edge of the high energy region is smeared by the 
missing invariant mass cut. It is interesting to note that, unless the branching ratio is not close 
to unity, we can find the two £™ ax edges with m G / 7^ and m G > = 0, which can determine 
both m x and m G i. Moreover, we can also determine the branching ratio from the shape of the 
distributions, i.e. the information on the coupling. 




(33) 
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Figure 4: The distributions of the leading photon energy (left) and of the missing invariant mass 
(right) for e + e~ — > XX ~ > 77 & t y/s = 500 GeV, where the pseudo-goldstino mass is fixed at 85 
GeV (red) and 125 GeV (blue). The mc = case (black) is also shown as the reference point as well 
as the SM background. 

3.4 Goldstini production in pp collisions 

We now turn to consider the processes which are relevant to the LHC. 

Similar to the process e + e~ — > \G' in Sec. I3.2[ the cross section of pp — > \G' is proportional 
to K 2 /m^ 2 and could be enhanced by the factor K. However, rather large K values of O(10 4 ~ 5 ) 
are needed to obtain a visible cross section for an eV-order gravitino, leading to a bound for 
the K values by the 7 +I^t events at the Tevatron [J8] similar to the LEP bound discussed 
above. We note that at the parton level the amplitudes are the same as the ones studied in the 
previous section after replacing the incoming electrons with the quarksfl 

The clean 77 signal is given by the neutralino LOSP pair production. However, the 
cross section is too small Ppp-^xx. ~ 

0.3 (1.2) fb at y/s = 7 (14) TeV at SPS8) to be significant 
over the SM background^ It is worth to emphasize again that the emitted photons associated 
with a pseudo-goldstino are softer than those with a true goldstino, which makes it difficult to 



4 Note that only in this process, and the following one, we would need to specify some squark masses. These 
can also be easily extracted from SPS8, however as we will argue the most interesting process that we study in 

more details does not involve intermediate coloured superpartners. 

5 It is well known that the NLO QCD corrections enhance the LO cross section by a factor of 1.3-1.4 for 
y/s = 14 TeV [49] . It is also interesting to note that gg collisions can give a certain contribution to the neutralino 
pair production through one loop [50j . The corresponding amplitudes are suppressed by the loop factor, while 
they are enhanced because of the larger gluon PDF, and also because all (s)quarks can run in the loop. 
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p T of the leading photon [GeV] p T of the leading photon [GeV] 

Figure 5: The pt distributions of the leading photon for pp — > — > l + l~ +77+$r at yfs = 14 TeV 
for rriQi = 85 GeV (left) and 125 GeV (right). The rag' = case is also shown as a reference with 
black lines. 

apply some cuts to enhance the signal over the background. In other words, the experimental 
constraints for the SPS8 point in the standard minimal GMSB model as well as for general 
gauge mediation [23] (e.g. [51j|52l[20j[53j|5l]) could be eased in the multiple goldstini scenario. 

Among the exclusive processes with at most two extra particles in the final state, the cleanest 
one is the one where the two photons and missing energy are accompanied by two leptons. By 
far the main contribution to the signal comes from the pair production of sleptons, which 
subsequently decay into a lepton and the neutralino LOSP, pp — >■ l^/ L l^/ L — > + 77 +$t 
(/ = e, fi). The electron and muon pairs give the same contributions due to the degeneracies 
between the first two slepton families. The Standard Model background is completely negligible 
compared to the signal of new physics. 

The signal cross section is 1.1 (3.3) fb at y/s = 7 (14) TeV at SPS8, where the slepton masses 
are m i R/L — 180.2/358.2 GeV and we employ the CTEQ6L1 PDFs [55] with the factorization 
scale chosen as \i = (mf + mi )/2^| Here, the leptons and photons are required to have 
Pt > 20 GeV, \rj\ < 2.5, and Ru m > 0.4, where pt and rj are the transverse momentum and 
the pseudorapidity of a final-state particle, respectively, and Rij describes the separation of the 
two particles in the plane of the pseudorapidity and the azimuthal angle. 

Because of the peculiarities of hadronic collisions, raising the center-of-mass energy from 

6 The NLO cross section is about 1.35 times larger than the LO one at y/s = 14 TeV [SHI HH]- The gluon 
fusion contribution to the slepton pair productions has been also studied in |571 158] . 
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Figure 6: The same as Fig. [5J but the missing transverse energy distribution. 



7 TeV to 14 TeV increases the total cross section but does not alter significantly the shape of 
the distributions. Therefore, we present the results only for the 14 TeV LHC here. In Figs. 
and [6J the distributions of the px of the leading photon and of the transverse missing energy 
are shown, respectively, where the same benchmark points are taken as in Sec. 13.31 The SM 
background is invisibly small. The two plotted variables show different shapes with respect 
to the standard goldstino scenario (mgi = 0) depending on the mass and the coupling of the 
pseudo-goldstino. This fact could in principle be used to extract the masses and the couplings 
by using techniques explained in [591 EQl EU [62]. We note that in the case of mc = 125 GeV 
the signal cross section is largely reduced by the experimental cuts, especially for large K^. 
This is due to the fact that, for more massive pseudo-goldstino, the emitted photons will be 
softer, and hence more excluded by the kinematical cuts. 

Before closing the section, we note that the production cross sections of the colored SUSY 
particles at SPS8 are small due to their large masses in the TeV range. In general, although we 
have presented the exclusive signals without jets here, the inclusive search is also interesting 
and will be reported elsewhere. 



4 The case of the stau LOSP 

We now move on to the case where the f is the LOSP, and use the Lagrangian f l22|) to study 
it. The effect of the coefficients K TL and K TR is that of enhancing the G' production. 
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Figure 7: Branching rations of f — > tG' at the SPS7 benchmark point, for various values of K TL and 



K 



tr ■ 



4.1 Decays of the stau 

The decay amplitude squared for the process f — > tG' can be easily computed in terms of the 
coefficients in the Lagrangian yielding 

,2N 2 



\M\' 



m% — ml 



{{K 2 TL cos 2 6 f + K 2 TR sin 2 9 f )(m 2 f - m 2 G , - m 2 T ) 
+ 4K TL K TR sin f cos 9 f m G /m T } . 



The width can be written as 



Y(f -)• tG') 







16iTmf 



(34) 



(35) 



with j3 = /3(^f, — y") as in ( )27|) . Notice that the width depends not only on the magnitude 

T T 

but also on the ratio of the coefficients K TL and K TR . The well known expression for the true 
goldstino case is obtained from the above by simply setting K TL = K TR = 1 and tuqi = 0. 

For illustration we plot the branching ratios of the decay f — >■ tG', choosing the stau 
mass and mixing angle given by the SPS7 benchmark point, namely m? = 124.0 GeV and 
cos 9f = 0.154. The smallness of cos#f at the SPS7 point makes the branching ratio largely 
independent from K TL but this would of course change for a different mixing angle. 



19 




500 



p T of the leading tau [GeV] Missing transverse energy [GeV] 

Figure 8: The distributions of pt of the leading tau (left) and of the missing transverse energy (right) 
for pp — > f + f~ — > t + t~ +]pT at -y/i = 14 TeV for mc 1 = 85 GeV. The mg/ = case (black) is also 
shown as the reference point as well as the SM background. 

The total width of the stau LOSP at SPS7, where ma(— r/13/2) = 0.76 eV, is of the order 
of 10 -10 GeV for K TL = K TR = 1, and hence the stau will decay promptly into the detector. 
Scenarios with long lived staus could also be envisioned, with even the possibility of stopping 
them and then measuring their decays. In that case the branching ratio and the mass of 
the invisible particle(s) should be more directly accessible. We will not delve further on this 
scenario. 

4.2 Goldstini and di-tau production 

We now turn to stau production at colliders and subsequent decay into taus and (pseudo)- 
goldstini. Though we could first consider stau pair production at e + e~ colliders, we will omit 
this rather straightforward exercise. In the neutralino LOSP scenario, the benefit of this set up 
was the ability to infer the neutralino and pseudo-goldstino masses from the energy distribution 
of the emitted photons. In the present case the same considerations are practically harder to 
achieve because of the difficulties inherent in tau reconstruction. We thus proceed to consider 
directly the case of hadron colliders. 

The process we are interested in is pp — > f + f~ — > t + t~ +$t- The lighter stau-pair 
production cross section at SPS7 is 41 fb at y/s = 14 TeV0 while the irreducible SM background 
t + t~vv, mainly coming from Z and W pair productions, is ten times larger with the minimal 

7 See footnote El and also [63] for more details on the stau-pair production. 
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rag' = 


tug 1 = 85 
K T = 1 K T = 2 K T = 10 


background 


a [fb] 
s/b 

S (100 (30) fb' 1 ) 


15.5 
0.18 
5.1 (2.8) 


13.2 10.1 5.9 
0.15 0.12 0.07 
4.4 (2.4) 3.4 (1.8) 2.0 (1.1) 


86.4 



Table 1: Cross sections for the signals and the background. The significance S in fl36l) is 
calculated with an integrated luminosity of 100 (30) fb -1 and a di-tau detection efficiency 
A = 0.1. 

cuts px T > 20 GeV, \rj T \ < 2.5, and R TT > 0.4. Since the background is located in the low tau- 
tau invariant mass region, we impose the additional cut M TT > 150 GeV in the analyses below. 
Figure [8] shows the distributions of the pt of the leading tau and of the missing transverse 
energy for the following benchmark points covering different branching ratios for stau decays 
into G or G": 

la. m G i = 85 GeV with K TL = K TR = 1 [B(f -> tG') ~ 0.2] 
lb. m G i = 85 GeV with K TL = K TR = 2 [B(f -> rG 1 ) ~ 0.5] 
lc. m G i = 85 GeV with K TL = K TR = 10 [B(f -> tG') ~ 1] 

The reference point mc = 0, i.e. the standard goldstino scenario, is also shown by a black solid 
line. As in the neutralino LOSP scenario, a sizable branching ratio of stau decays to pseudo- 
goldstini has an impact on the shape of the distribution, once again making the spectrum 
softer. Another noticeable fact is that the new physics signal becomes dominant especially for 
high missing transverse energy only when the stau decays to true goldstinos are significant. 
Therefore, the pseudo-goldstini scenario makes it harder to achieve enough significance for its 
observation. 

In order to quantify the sensitivity of the missing transverse energy cut, in Tables [1] and [2] we 
list the cross sections, signal (s) over background (b) and significances (S) for our benchmarks, 
without and with a cut on If)? at 100 GeV respectively. Taking into account the branching 
ratio of the hadronic tau decays, B{Th a d) = 0.648 [2Z|, and the efficiency of the hadronic tau 
identification, e = 0.5 [6j|,[65], we assume the di-tau detection efficiency A = (B(rh a d) x e) 2 ~ 
0.1. We use the signal significance defined as [66] 

S = ^2((s + b) \n{l + s/b) -s). (36) 
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rag' = 


tug 1 = 85 
K T = 1 K T = 2 K T = 10 


background 


a [fb] 
s/b 

S (100 (30) fb' 1 ) 


5.6 
1.59 
7.9 (4.3) 


4.8 3.4 0.9 
1.36 0.95 0.25 
6.9 (3.8) 5.0 (2.7) 1.4 (0.8) 


3.5 



Table 2: The same as Tabled], but with the additional $t > 100 GeV cut. 



We observe that for 100 fb -1 of integrated luminosity, a satisfactory significance can be achieved 
for most of our benchmarks, especially for Ifi? > 100 GeV. However, in the K T = 10 case, i.e. 
in the case of the heavy pseudo-goldstino with the enhanced coupling, the cut on ]px does not 
improve the significance of the signal. 

Finally we note that, although we took K TL = K TR for simplicity, tau polarization may be 
exploited to determine the case that the two coefficients are different, K TL ^ K TR [67] . 

5 Conclusions 

In this paper we have analyzed how the traditional expectations from a low scale gauge mediated 
model can be modified by the presence of a massive pseudo-goldstino NLSP in a scenario with 
multiple-sector SUSY breaking. 

We first considered the case of a gaugino-like neutralino LOSP. We showed that the decay 
modes of the LOSP into a photon or Z-boson and a pseudo-goldstino can be significant. We 
studied in details the goldstini phenomenology in the photon(s) plus missing energy signals in 
e + e~ and pp collisions. Our aim was to provide clues to interpret prompt photon plus missing 
energy signals at the LHC We found that the resulting photon spectrum is typically softer 
and with different shapes, compared to the standard gauge mediation scenario with only one 
hidden sector. 

We then proceeded to consider the case of a stau LOSP. Similarly to the previous case, 
the stau can significantly decay into a tau and a pseudo-goldstino. We studied a possible tau 
pair plus missing energy signal at the LHC, again finding that the pseudo-goldstino scenario 
leads to an altered shape in the tau energy distributions. The signature however seems more 
promising in the neutralino LOSP scenario, not only because photons have cleaner experimental 
features than taus, but also because the SM background is virtually absent in the process that 
we studied. 
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Although we performed our analyses at the particular benchmark points (SPS8 for \ LOSP 
and SPS7 for f LOSP) to make a comparison with the standard gauge mediation model, our 
goldstini model can be easily extended to other benchmark points and LOSP scenarios due 
to the implementation into a event generator, MadGraph 5. We emphasized the dependence 
of the experimental signatures on the new parameters, namely the pseudo-goldstino mass and 
its couplings. A rather model independent feature of the pseudo-goldstini scenario is that it 
provides a scenario which is fully consistent with gauge mediation of supersymmetry breaking, 
while it softens significantly the energy of the SM products of the LOSP decays, due to the 
conspicuous mass of the pseudo-goldstino. In this way such a scenario evades most experimental 
bounds that tend to exclude gauge mediated scenarios. 

We would like to suggest that models like the one we have presented here, where the missing 
energy is carried away by two different particles, one rather massive and the other almost 
massless, and with similar couplings but with different strength, deserve to be considered in 
LHC data analysis. 
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